Uncertainties in the Anti-neutrino Production at Nuclear Reactors 
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Anti-neutrino emission rates from nuclear reactors are determined from thermal power measure- 
ments and fission rate calculations. The uncertainties in these quantities for commercial power 
plants and their impact on the calculated interaction rates in v e detectors is examined. We discuss 
reactor-to-reactor correlations between the leading uncertainties, and their relevance to reactor v e 
experiments. 
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I. INTRODUCTION 



Electron anti-neutrinos from large commercial nuclear 
reactors are playing an important role in the exploration 
of neutrino oscillations The choice of distance be- 
tween source and detector allows one to conveniently tune 
an experiment's sensitivity to either the atmospheric or 
solar neutrino mass splitting. KamLAND 0, H|> 0] has 
established a connection between the MSW effect in the 
sun [5[ and vacuum oscillations with anti-neutrinos, and 
has provided the most accurate measurement of Am^ 
to date. Earlier, CHOOZ [| and Palo Verde 0] pro- 
vided what are still the best upper limits on the mixing 
angle # 13 . A variety of new experiments are being con- 
structed to perform more sensitive measurements of the 
mixing angle 613 [1, [l(| [HI [HJ . Anti-neutrino detec- 
tors at very short distance are being explored with the 
purpose of detecting coherent neutrino-nucleon scatter- 
ing (CiaNS) [TH, and for testing plutonium diversion at 
commercial reactors [l4| . A proper understanding of the 
systematic effects involved in the modeling of the reactor 
anti-neutrino flux and energy spectrum is essential for 
these experiments. Moreover, as most experiments mea- 
sure D e s from more than one reactor, it is important to 
understand not only the magnitude of these effects, but 
also the correlations between uncertainties from different 
sources. 

Except for the case of C^NS, which will not be ad- 
dressed further here, reactor anti-neutrinos are usually 
detected using the inverse-/? reaction v e + p — > n + e + 
whose correlated signature helps reduce backgrounds but 
is only sensitive to v e s with energy above 1.8 MeV. The 
kinetic energy of the positron gives a measure of the in- 
coming v e energy. Neutrino oscillation parameters may 
be extracted from the data by fitting the observed v e 
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spectrum j^- to the following equation: 
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Here, N p is the number of target protons, t{E v ) is the 
energy-dependent detection efficiency, a(E u ) is the detec- 
tion cross-section, P ee (E u , L^) is the oscillation survival 
probability for P e s traveling a distance Lfc from reactor 
k to the detector, and Sk{E u ) is the v e spectrum emit- 
ted by reactor k. The time variation of the flux may 
be included in the fit to discriminate the time-varying 
reactor signal from constant backgrounds. Equation [1] 
should technically be multiplied by the detector resolu- 
tion function and integrated over E v , but this detail is 
inconsequential for our analysis. 

Nuclear power reactors operate on the principle that 
the fission of U and Pu isotopes and the subsequent de- 
cays of their daughter fragments release energy, gener- 
ating heat. Large Q- value /3-decays of unstable fission 
fragments are primarily responsible for the v e emission 
of nuclear reactors. Decays of long-lived isotopes in the 
nuclear fuel and in spent fuel elements stored at the reac- 
tor site contribute at the sub-percent level to the detected 
D e flux; this contribution has been treated elsewhere (l5j 
and is beyond the scope of this paper. Considering only 
fission reactions, the final term in Equation [T] may be 
expanded as 
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where dN vi jdE v is the v e emission spectrum per fission of 
isotope i, and fi is the number of fissions of isotope i dur- 
ing the data taking period. Since >99.9% of the energy- 
and v f .-producing fissions are from 235 U, 238 U, 239 Pu, and 
241 Pu [16| (see Figure [lj, the summation in Equation [2] 
is performed over just these 4 isotopes. For 235 U, 239 Pu, 
and 241 Pu, the emitted j? e spectra dN v JdE v are derived 
from /3-spectrum measurements of the fissioning of the 
isotopes by thermal neutrons [lj], [3 • The uncertainty 
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in the v & spectral shape derived from these measurements 
was investigated recently in [Iff] . For 238 U, no measure- 
ments are available, so theoretical calculations of its v e 
emission must be used (2pj . 

The fi can be obtained from detailed simulations of 
the reactor core throughout the data taking period. The 
output fi must obey the thermal energy constraint 
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where Wth is the total thermal energy produced by the 
reactor during the time period considered, and denotes 
the energy released per fission of isotope i. To reduce sen- 
sitivity to errors in the simulation codes, the codes are 
typically used to obtain not the fi directly but the fis- 
sion fractions fi/F, where F = Y2%fi- Measurements of 
the total generated thermal power taken regularly during 
reactor operation may then be used to obtain the total 
numbers of fissions of each isotope using Equation [3] 

The appropriate simulation codes are cumbersome to 
run, and are often proprietary, and therefore are not al- 
ways available for direct use by scientific collaborations. 
Obtaining the fi/F therefore requires special agreements 
with reactor operators. Depending on the terms of the 
agreement, fine grained time-dependent fission fractions 
may be obtained. Fine-grained data for Wth are com- 
monly available, since such information must be made 
available to nuclear regulatory bodies. Reactor v e ex- 
periments then use this information to calculate the P e 
signal using an equation of the following form: 
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where for simplicity we have assumed a single reactor, 
and have integrated over F v . The terms in Equations [T] 
and [2] not appearing in Equation |4] have been grouped 
into the rij, the number of detected v e per fission of iso- 
tope i. 

Most recent reactor v e experiments, including Kam- 
LAND i, i, ||, CHOOZ @, Palo Verde H and 
Bugey [21], have used Equation^ or variations on it, to 
perform their signal calculations. In the first two exper- 
iments, the systematic uncertainty in the signal calcula- 
tion was taken to be the quadratic sum of the uncertainty 
in the reactor power measurements which yield Wth and 
some estimate of the error due to the uncertainty in the 
fission fractions fi/F. In the case of Palo Verde this lat- 
ter component was estimated from comparisons between 
calculated and measured isotopic concentrations in spent 
fuel elements [l6[ . In the case of KamLAND this compo- 
nent was taken from comparisons of a simplified reactor 
model to detailed simulations 0, H, d, , which assumes 
that the uncertainty of the simplified reactor model sig- 
nificantly exceeds that of the detailed simulations. How- 
ever, no systematic uncertainty was assigned to the de- 
tailed reactor simulations. The CHOOZ experiment 0| 



leveraged the short baseline Bugey observation to limit 
the reaction cross section uncertainty, scaling their re- 
sult to agree with the Bugey-measured cross section per 
fission. This treatment allowed CHOOZ to significantly 
reduce the error associated with the v e emission spectra 
(dN v /dE v ). However, their method is legitimate only 
to the extent that the reactor simulations employed by 
the CHOOZ and Bugey reactor operators reliably model 
the different reactors with their specific fuel compositions 
and operation histories. In all cases, the combined sys- 
tematic error was estimated to be within 2-3%. At this 
level, the different methods employed by each experiment 
are likely not overly aggressive. 

With detector-based uncertainties in current and next- 
generation reactor v e experiments approaching the ~1% 
scale, a more careful treatment of these reactor-specific 
uncertainties is necessary, particularly for experiments 
with multiple reactor sources in which correlations be- 
come important. In this paper we outline such a detailed 
treatment of these uncertainties and demonstrate its ap- 
plication to a counting analysis of D e 's from single and 
multiple commercial reactor sources. In Section|TT]we dis- 
cuss in detail the uncertainty in reactor thermal power 
measurements. In Section IIIII we examine the Monte- 
Carlo estimation of the fi and, expanding on the treat- 
ment in [l6j and [23| , derive uncertainties for these based 
on a large body of spent fuel isotopic concentration com- 
parisons used to verify the codes. In Section HVl we com- 
bine the uncertainties from Wth and the fi to estimate 
their contribution to the uncertainty in the anti-neutrino 
yield. We also discuss the applications of our calculation 
to specific experimental configurations, including multi- 
ple reactor sources. We draw our conclusions in Sec- 
tion El 

Our analysis addresses issues common to large com- 
mercial power reactors, especially Pressurized- Water- 
Reactors (PWRs) and Boiling- Water- Reactors (BWRs), 
using low-enrichment fuel. Explanation of the operation 
of PWRs and BWRs may be found in [24[ . Modern neu- 
trino experiments are almost exclusively using v e from 
such reactors. Minor contributions from other reactor 
types are not considered here. All uncertainties are given 
at the 68.3% confidence level. 



II. THERMAL POWER UNCERTAINTY 

The most accurate measurement of a reactor's ther- 
mal power is given by a calculation of the energy bal- 
ance around the reactor vessel (BWR) or steam genera- 
tor (PWR). This requires accurate measurements of feed- 
water mass flow and temperature, steam enthalpy and 
moisture content, and reactor coolant-cycle heat gains 
and losses 25]. For the PWRs at Diablo Canyon Power 
Plant (DCPP) (California, USA), the thermal balance is 
written as 

Qc — Qs + Qltnd + Qr&c — Qrcp — Qpz- (5) 
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Qc is the core thermal output, the time- integral of which 
gives W t h- The largest component of it, Q s , is the 
power extracted from the steam produced directly in the 
steam generator. Smaller corrections are represented by: 
Qltnd, the power lost in the water clean-up system; 
QiUzCi the power losses to the external environment due 
to radiation and convection; Que Pi the contribution to 
Qs from heating of the working fluid by the circulation 
pumps; and Qpz, the contribution to Qs from heating 
of the working fluid by the pressurizers. The correction 
terms together account for 0.3-0.4% of the total Q c . 

A similar thermal balance may be written for other 
PWRs as well as for BWRs. In the case of BWRs, the 
main component is also in the form of steam produced 
by the reactor. For both reactor types, Qs is evaluated 
according to 

Qs = m s (h out - h in ), (6) 

where m s is the mass flow rate of the feed-water to the 
steam generator (PWR) or reactor vessel (BWR), and 
h ut — hi n is the specific enthalpy rise in the steam gen- 
eration. We will discuss the uncertainties in these terms 
below. The uncertainties on the terms other than Qs 
in the thermal balance equation (Equation [5]) contribute 
negligibly to the uncertainty of the core thermal power 
and will not be addressed further here. 

The enthalpy rise is calculated, using steam tables i 26] , 
from inlet and outlet values of the pressure and tempera- 
ture, and from the moisture content on the secondary side 
of the steam generator or reactor vessel. The enthalpy 
uncertainty has several contributions deriving from mea- 
sured and calculated quantities. The errors in the tem- 
perature and pressure measurements are assumed to be 
random after correction for known systematic contribu- 
tions. The moisture content of the steam is generally 
known to grow during the operation cycle because of re- 
duced moisture removal efficiency of the steam separa- 
tors due to a slow deposition of eroded metal particles. 
The moisture content uncertainty of the steam may be 
treated as a systematic uncertainty since all reactors un- 
dergo this aging process. The calculation of the enthalpy 
rise from these inputs using steam tables contributes an 
additional systematic uncertainty of <0.2% [27l |. 

As an example, Table Q] shows the uncertainties of the 
quantities relevant to the enthalpy calculation for DCPP. 
The moisture content at DCPP's steam generators grows 
with time, as mentioned above; the number in Table [I] is 
an average over the period of operation. Over longer pe- 
riods the moisture content may grow, as explained above, 
by a factor of two over the level represented in Table [H 
but even then it will comprise a sub-dominant contribu- 
tion to the error budget. A similar analysis for Beaver 
Valley Unit 2 yielded enthalpy-related contributions to 
the thermal power uncertainty totaling 0.16% (28|. A 
more general discussion in [29( uses an enthalpy uncer- 
tainty of 0.24%. For a generic reactor or when a more de- 
tailed analysis is unavailable, we suggest conservatively 
assuming 0.15% random uncertainty (for the p and T 



TABLE I: Uncertainty contribution to the enthalpy from mea- 
sured input quantities at DCPP. The change of the enthalpy 
vs given quantityUAEnthalpy/AQuantity), is calculated us- 
ing steam tables [26j . pt n , Ti„ and p ou t are the inlet pressure, 
the inlet temperature, and the outlet pressure, respectively. 
The present moisture carryover value at DCPP is assumed in 
the table. The quadratic sum conservatively assumes a full 
0.2% error related to the use of steam tables to calculate the 
enthalpy rise. 



Typical Quantity AEnthalpy Enthalpy 
Value Error [%] AQuantity Error [%] 



Pin 


6.9 MPa 


0.50 


0.002 


0.001 


Ti n 


221°C 


0.12 


1.153 


0.138 


Pout 


5.6 MPa 


0.94 


0.018 


0.017 


Moisture content 


0.99 


0.05 


0.562 


0.028 


Steam tables 








<0.2 


Quadratic Sum 








0.25 



measurements) and 0.2% systematic uncertainty (due to 
the moisture content and use of steam tables) for the 
enthalpy rise. 

The uncertainty in m s gives the largest contribution to 
the error in the thermal power [30L HH . In order to pro- 
vide as accurate a measurement as possible, the flow rate 
is measured in a long, straight section of pipe where fully 
developed turbulent flow is established. Different types 
of flow meters may be employed there. Traditionally the 
system is instrumented with nozzle-based flow meters op- 
crating on the principle of the Venturi effect. An early 
analysis yielded a total flow uncertainty for these me- 
ters of 1.7%, of which 1.6% is attributable to random 
errors, while the systematic uncertainty is determined 
to be 0.5% [Hj]. According to more recent experience 
at DCPP, properly calibrated and maintained Venturi 
meters can be operated with an initially low total un- 
certainty of about 0.7%. However, the accuracy of the 
Venturi meters fluctuates with time due to material de- 
position, or "fouling" [3(3], which leads to smaller cross 
sectional area for the flow through the device and results 
in indicated flow rates exceeding the true value. Depend- 
ing upon the nozzle style and the pH factor of the feed- 
water, this fouling can grow as high as 3% over periods 
of a few months to a few years. In-situ comparisons be- 
tween Venturi meters and more modern instrumentation 
(discussed below) shows an average systematic difference 
of +0.6%, and a statistical spread of -1.4% [H,!!]. 

Fouling in Venturi flow meters results in economic 
losses for the nuclear power industry. Nuclear reactor 
regulations in the US and elsewhere require operators to 
apply a safety margin accounting for the uncertainty in 
the thermal power measurement when performing loss- 
of-coolant accident and emergency core cooling system 
performance analyses. However, the regulatory practice 
varies from country to country [35j . The size of the 
safety margin is 2% in the US, but may be reduced if the 
thermal power measurement uncertainties can be demon- 
strated to be smaller than that level [36[. For reactors 
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instrumented with Venturi meters, for which lower un- 
certainties cannot be demonstrated as a result of fouling, 
this in effect leads to operation at 2% below the plant's li- 
censed operating power limit. A similar 2% safety margin 
is required in Japan [371 ] . Moreover, fouling itself results 
in an overestimation of the thermal power, leading to ac- 
tual operation up to a few percent below the indicated 
level, and thus a lower electric power output. This two- 
fold impact of fouling gives the nuclear power industry a 
high incentive to deploy a replacement technology with 
smaller, well-understood uncertainties. 

Using a new generation of ultrasonic flow meters 
(UFM) not affected by fouling, the thermal power uncer- 
tainty may be significantly reduced. UFMs use an elec- 
tronic transducer with no differential pressure elements. 
There are currently two types of UFM systems routinely 
used in the nuclear electric generation industry to accu- 
rately calculate feed- water flow in high purity water [38| . 
One style is the transit-time meter, while the other is 
the cross-flow meter. The transit-time UFM emits ultra- 
sound signals diagonally through the fluid upstream and 
downstream along the same path. The measured time 
difference between the two paths is proportional to the 
velocity of the fluid in the pipe. Accuracies of 0.2-0.5% 
are reported [H, H(| for transit-time type UFMs. The 
cross-flow UFM measures the time taken by a unique 
pattern of turbulent eddies in the fluid to travel between 
two pairs of transducers at some known distance apart 
along the pipe. The cross-flow devices are mounted exter- 
nally to the pipe. The cross-flow UFM with its associated 
hardware and software is able to achieve a flow measure- 
ment uncertainty of 0.25% or better, although the ability 
to achieve these uncertainties in real reactor conditions 
has been questioned [4l[ . This type of meter is employed 
at DCPP, with an uncertainty evaluated in-situ within 
the 0.4-0.7% range. The principles of operation of both 
transit-time and cross-flow UFMs are described in detail 

in m. 

The measurement uncertainties of transit-time UFMs 
and their breakdown into statistical and systematic com- 
ponents were studied in detail by Estrada [4(| of Caldon 
Inc., a flow meter manufacturer for the nuclear industry. 
This work found total uncertainties in the measured mass 
flow rates of 0.45% for models with externally mounted 
transducers, and 0.20% for intrusive-type models with 4 
pairs of transducers. We refer to these as Type I and 
Type II UFMs, respectively. These uncertainties origi- 
nate from the differences between the water axial veloc- 
ity profile in the test facility and in the measurement 
system in a plant, the uncertainty in the measurement 
of the ultrasound acoustic path in the water pipe, the 
imperfect knowledge of the dimensions of the measure- 
ment systems, and the uncertainty in the measurement 
of the time of flight of the acoustic pulses, including non- 
fluid delays. The water flow uncertainties are divided 
into random and systematic contributions in Table [II] 

We note that the 0.45% uncertainty based on Estrada's 
study is consistent with the 0.4-0.7% uncertainty esti- 



TABLE II: Typical flow meter uncertainties. Type I corre- 
sponds to transit-time (TT) UFMs with ultrasonic transduc- 
ers externally mounted to the feed-water pipe. Type II cor- 
responds to transit-time UFMs with ultrasonic transducers 
in-line with the pipe via a spool piece that integrates four 
pairs of ultrasonic transducers, forming four chordal paths. 
Type I-II UFM water flow uncertainties are based on Caldon 
Inc. experience [4(j. The flow meter employed at DCPP is of 
the cross-flow (CF) type. 



Flow Meter 


Random [%] Syst. [%] 


Total [%] 


Venturi 


1.4 0.6 


1.5 


Type I TT UFM 


0.2 0.4 


0.45 


Type II TT UFM 


0.09 0.18 


0.20 


DCPP CF UFM 


0.3-0.6 0.1-0.2 


0.4-0.7 


TABLE III: Typical thermal balance quantities for the Diablo 


Canyon PWR reactors. 




Quantity 


Typical Value 


Error [%] 




1887 kg/s 


0.4-0.7 


A/l = h ou t ~ hin 


1.819 ■ 10 6 J/kg 


0.25 


Qs 


3433 MWt 


^0.7 


Qltnd 


1.83 MWt 


nil 


Qr&c 


0.65 MWt 


nil 


Qrcp 


14 MWt 


nil 


Qpz 


0.21 MWt 


nil 


Qc 


3421 MWt 


^0.7 



mated by DCPP for cross-flow meters. Several other as- 
sessments of mass flow rate uncertainties, also consistent 
with these values, may be found in the literature. The 
generic iscussion in [29j, for example, quotes feed- water 
flow uncertainties of 0.4%. For a generic UFM, we will 
use the values listed for the Type I transit-time UFM, 
which has the largest estimated systematic uncertainty, 
and gives a total error consistent with values considered 
by the US Nuclear Regulatory Commission in a discus- 
sion on the use of UFMs [4l|. 

The errors on the enthalpy rise and m s may be added 
in quadrature to obtain the full error on Qs- Combin- 
ing the mass flow uncertainty (as given in Table [ill for 
Type I transit-time UFMs) with an enthalpy uncertainty 
of 0.25%, we get a total thermal power uncertainty of 
0.51%, of which 0.45% is systematic and 0.25% will vary 
from reactor-to-reactor. Typical values for the quanti- 
ties describing the heat balance and their error are given 
in Table UTT1 for DCPP, using two PRWs equipped with 
cross-flow type UFMs. An internal DCPP study char- 
acterized about a quarter of the thermal power uncer- 
tainty as correlated. The error in the determination of 
Qs, which is the only non- negligible component of Qc 
and hence determines the uncertainty on Wth, is domi- 
nated by the uncertainty in m s . For precise evaluation of 
the errors, the enthalpy uncertainty has to be taken into 
account. 

It is worth mentioning that most feed-water flow ap- 
plications require two or more flow measurements. In 
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PWRs, for example, the flow is measured in each steam 
generator. In BWRs, the flow is normally measured 
in each of two main feed headers. In these cases, the 
random uncertainties may be reduced by a factor of 
1/V num ber of generators or l/y/2, respectively. 

It is also worth noting that estimates of flow rate mea- 
surement uncertainties to date have been performed in 
laboratory settings at Reynolds numbers of up to 10 6 . 
However, Reynolds numbers in actual power plants are 
as high as 10 7 [5^, [43|, |44[ ■ The uncertainty estimates de- 
scribed above are therefore extrapolations from the lower 
Reynolds number evaluations. Although the error asso- 
ciated with these extrapolations is expected to be small, 
new test facilities have been proposed to generate real- 
istic flow conditions with a higher Reynolds number to 
test these extrapolations [451 ]. 



III. UNCERTAINTIES IN FISSION 
CALCULATIONS 
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During the power cycle of a nuclear reactor, the com- 
position of the fuel changes as Pu isotopes are bred and 
U is depleted. At the end of the power cycle, some frac- 
tion of the fuel is replaced, and the remaining fuel el- 
ements are shuffled to optimize burnup. Knowledge of 
the detailed, time-dependent radionuclide content of the 
spent fuel throughout the power cycle is of interest to 
reactor operators and regulatory groups because it im- 
pacts shielding requirements, dose rate analysis, toxicity, 
waste storage and handling considerations, and critical- 
ity safety [4y|. To obtain this knowledge, reactor opera- 
tors employ Monte-Carlo simulations of the reactor core. 
The simulations begin with an initial fuel composition 
and a physical model of the core assembly. An initial 
neutron flux solution for the core is computed based on a 
number of inputs, including the measured thermal power 
generation, and other operating parameters such as the 
pressure, temperature, and flow rate of the cooling sys- 
tem, neutron moderation parameters, etc. Cross-section 
and decay data libraries are used to compute the rate of 
various interactions and update the fuel composition as 
a function of position. This process is iterated in a series 
of time-steps throughout the reactor power cycle. 

As a by-product, these reactor core simulations may 
be made to output the number of fissions fi needed for 
the estimation of the signal in reactor anti-neutrino ex- 
periments. A plot of the time-dependent fission rates 
calculated for a typical power cycle of one of the reactors 
at the Palo Verde Nuclear Generating Station (PVNGS) 
in Arizona, USA, is shown in Figure [TJ Integrating over 
the power cycle gives the fi required by Equation [4] to 
compute the v e signal. 

Uncertainties in these simulation codes originate from 
a variety of sources [47j . These sources include the un- 
certainties in the input parameters, uncertainties in the 
nuclear cross-section and decay data used by the codes, 
approximations made in the modeling of the reactor core, 



FIG. 1: Fission rates for various isotopes during a typical 500 
day fuel cycle for one of the Palo Verde Nuclear Generating 
Station reactors in Arizona (USA). >99.9% of the z? e s are 
emitted by 235 U, 238 U, 239 Pu, and 241 Pu 0. 



and numerical approximations used in the computations 
themselves. Reference [l6j explored the input parameter 
uncertainty contribution with the reactor core simula- 
tion code used at PVNGS [H, HU by calculating the v e 
signal variation for a given deviation of several of the 
input parameters. For the PVNGS simulation, the nor- 
malization to the measured thermal power represented 
by Equation 0] was done within the code. The code was 
run for a set of 500-day power cycles with one parame- 
ter varied at a time. The results of the study are shown 
in Figure H The thermal power measurement has the 
largest effect, with a slope of 0.96. 

The most common means reactor operators use to 
globally assess the performance of the codes is to com- 
pare measurements of the isotopic content of spent fuel 
elements with the values predicted by the simulations. 
Such verifications have been performed for a number of 
codes at a variet y of PWRs and BWRs in the US, Japan 
and Europe [50l l5ll. l52l l53l]. Samples are obtained from 
representative locations throughout the core, commonly 
selected from regions that are difficult to simulate ac- 
curately, such as boundaries and hot-spots. The assay 
technique employed is typically isotopic dilution mass 
spectroscopy, for which errors of 1-3% are typical for the 
isotopes of interest in more recent measurements, 
errors below 0.3% have been acheived [46| . 

Figure [3] shows an example comparison between mea- 
sured and calculated concentrations, in kilogram per met- 
ric tonne uranium (kg/MTU), of 235 U, 238 U, 239 Pu and 
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FIG. 3: Measured vs. calculated concentrations of 235 U, 238 U, 
239 Pu and 241 Pu as a function of burnup for fuel elements 
taken from the Takahama Unit 3 PWR [461 ]. See the text for 
a definition of the units. 



FIG. 2: The change in the expected i/ e -signal as a function 
of five inputs to the core simulation, taken from [l6l |. The 
numbers in the key are slopes of the fitted lines. The lower 
panel zooms in on the non-power inputs so that their variation 
can be observed. 



2 Pu as a function of burnup for a set of fuel elements 
processed at a typical PWR [46| • The quantity "burnup" 
in the abscissa is defined as the amount of energy (in 
GigaWatt days) extracted from a fuel element per unit 
initial mass of uranium (in MTU). The burnup is mea- 
sured to within ~3% for each sample using the 148 Nd 
method [55j . in which the integral number of fissions, de- 
termined from the sum of the 1 Nd fission product yields 
for the four isotopes of interest, is multiplied by the av- 
erage energy released per fission, weighted according to 
the calculated fission fractions. During a typical reactor 
cycle, fuel elements near the center of the reactor core re- 
ceive a higher burnup than those at the edges, due to the 
higher neutron flux at the center. This spatial variation, 
along with the number of cycles over which the samples 
were processed, is responsible for the different burnup 
values achieved in the samples plotted in Figure [3] Dur- 
ing normal operations, fuel elements are processed over 
several power cycles, and arc shuffled between each power 
cycle, until optimal burnup is reached for all fuel ele- 
ments. In practice, burnups up to ~50 GWd/MTU are 
achieved in the typical operation of modern commercial 
reactors [56| . 

Fractional differences in the heavy isotope concentra- 
tions Sci/ci were calculated from the measured (M) and 
calculated (C) concentrations using the following sign 



convention: 

5c _ c c - c M ^ 
c c M 

In Figure [4] we plot the values of 5a/ a (i — 235 U, 238 U, 
239 Pu, 241 Pu), using 159 comparisons of fuel element 
samples taken from ten PWRs and BWRs, modeled by 
a variety of core simulation codes. Details of the codes 
and reactors included in our analysis are listed with their 
references in Table IIVI Gaussian fits to the distributions 
are drawn for reference. The individual distributions for 
PWRs and BWRs (not shown) are equivalent within the 
available statistics. The concentration comparisons also 
do not exhibit a strong trend with burnup, as shown in 
Figure [5] 

Unfortunately, these isotopic concentration compar- 
isons only give indirect information on the uncertainty in 
the number of fissions. However, we can estimate 6fi/fi 
from Sci/ci by considering the nuclear processes responsi- 
ble for the changes in isotopic concentrations through the 
lifetime of a typical fuel clement. The U isotopes are the 
most straightforward to diagnose, since U breeding pro- 
cesses are negligible compared to U depletion processes. 
The decrease of the U concentrations with burnup results 
from a combination of fission and isotopic conversions due 
to other nuclear processes (e.g. n capture). The error on 
the total loss of U is the fractional error in the change in 
concentration after burnup B. We conservatively assign 
the error on the total loss of U entirely to fission pro- 
cesses, so that the fractional error on the number of U 
fissions is at most 

5_h 5(c v (0) - c v (B)) 

/u cv(Q)-c v (B) ' 1 ' 

where cu(0) and c\j(B) are, respectively, the initial and 
final concentrations of either U isotope after burnup B. 
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FIG. 4: The values Sc/c for 
calculated for each of the 159 comparisons of fuel element 
samples taken from a number of PWRs and BWRs modeled 
by a variety of core simulation codes 46 : ., 57*, 58, 59, :6(| . 
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FIG. 5: The values 8c/c for 235 U, 238 U, 239 Pu and 241 Pu as a 
function of burnup, calculated for each of the 159 fuel sample 
comparisons. 



The initial 235 U and 238 U concentrations (ctj(O)) are de- 
termined by the enrichment level of the fuel elements un- 
der consideration and are known very precisely. Hence 
the uncertainty <5(ctj(0) — c\j(B)) is dominated by the 
uncertainty in ctj(-B), allowing us to approximate 



Sfv 
fu 



1 



ScyjB) 
Cu (0)/ Cu (B)-iy cu(B) ' 



(9) 



where the last term is the 5c/ c plotted in Figured] As can 
be seen in Figure [3J which is typical for fuel processed in 
commercial PWRs and BWRs, the factor in parenthesis 
in Equation [5] is significantly different from one for 238 U 
but not for 235 U, except at low burnup. 

For the Pu isotopes the relationship between the 



TABLE IV: Codes and reactors included in our analysis. 
In cases where multiple tests were performed on the same 
code/reactor using older cross section files or neglecting some 
reactor configuration or operation details, we used the com- 
parisons for the most up-to-date data libraries and most de- 
tailed reactor model. 



rCCI. 


Reactor (Typo) W th [MWt] 




[46] 


Takahama Unit 3 (PWR) 


2652 


SAS2H 
HELIOS 


IP' 1 


Calvert Cliffs Unit l(PWR) 
H.B. Robinson Unit 2 (PWR) 


2560 
2192 




[58J 


Cooper (BWR) 
Gundrcmmingen (BWR) 
JDPR (BWR) 


2381 
801 
45 


SAS2H 


[59] 


Ti'ino Verccllcsc (PWR) 
Turkey Point Unit 3 (PWR) 


825 
2300 


SAS2H 


[60J 


San Onofre (PWR) 


1347 


SAS2H 


[61] 


Takahama Unit 3 (PWR) 
Fukushima Unit 2 (BWR) 


2652 
3293 


SWAT 
ORIGEN2 

SWAT 
ORIGEN2 


[62] 


Calvert Cliffs Unit l(PWR) 


2560 


ORNL-Assm 


TABLE V: Means (u) and standard deviations (<r) of Sfi/fi 


for commercial reactors. 






Fuel 


A* [%] 




a [%] 


235 tj 


0.09 ± 0.16 




2.0 


238 tj 


-0.81 ± 0.38 




4.8 


239 Pu 


0.74 ± 0.45 




5.7 


241 p u 


-0.32 ± 0.48 




6.0 



<5/pu//pu and the <5cp u /cp u is complicated by the pres- 
ence of significant breeding processes. However, the 
sources of uncertainty in <5/p u //pu are the same as those 
in <5cp u /cp u , indicating a similar magnitude. Moreover, 
since the Pu concentrations in fresh fuel used at typical 
commercial reactors starts at "zero" and the fission rate 
is at all values of B proportional to the concentration 
cpu(-B), we simply assume that <5/p u //pu ~ Scp u /cp u di- 
rectly. 

The four values of Sfi/fi (i = 235 U, 238 U, 239 Pu, 241 Pu) 
were calculated for each of the 159 comparisons of fuel 
element samples from the reactors and core simulation 
codes listed in Table IIV1 using Equation [5] for the U iso- 
topes and letting 5fp u /fp u sa i5cp u /cp u for the Pu iso- 
topes. The means and standard deviations of the dis- 
tributions for each isotope were extracted from Gaussian 
fits. The results are listed in Table fVl Note that with our 
sign convention defined in Equation [3 a positive Sfi/fi 
corresponds to an over-estimated fission count by the re- 
actor core simulations. 

Correlations between the Sfi/fi were studied as well. 
Figure [5] plots the Sfi/fi for pairs of fuel isotopes for 
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maintain the energy balance, the over-fissioning of this 
isotope must be accompanied by an under-fissioning of 
the other isotopes (and vice- versa). The fact that the 
largest anti-correlation is with 239 Pu, the next largest 
energy producing heavy isotope after 235 U, strengthens 
this argument. 

We also find a strong correlation between the two Pu 
isotopes, which can be explained by the fact that 239 Pu is 
a precursor to 241 Pu in one of three main isotope trans- 
formation chains within the reactor core [63| : 

238TJ ^239p u ^240p u ^241p u _»3«Cm 



(11) 



Since there is little or no Pu in a fuel element at zero bur- 
nup, over-fissioning of a Pu isotope is associated with its 
overproduction in the simulation via this isotope trans- 
formation chain. If one Pu isotope is overproduced, it is 
likely that the other will also be overproduced, resulting, 
in more fissions from both isotopes. 

The correlation between 238 U and Pu isotopes is a bit 
harder to diagnose. However, 238 U is only a minor player 
in the thermal power production, comprising only ~10% 
of the output, so a strong anti-correlation is not expected 
as it is for 235 U. And while 238 U is the parent of the 
transformation chain represented by Equation 111! that 
chain represents a series of neutron captures, whereas the 
correlations listed in Table PVTl correspond to fissions. The 
over-fissioning of 238 U does not a-priori imply a greater 
neutron capture rate to produce more of the Pu isotopes 
(so that more are available to fission), and vice- versa. 
The true 238 U-Pu correlation is a combination of these 
competing effects. 



FIG. 



Hi 



dependencies for 235 U, 238 U, 239 



Pu, and 241 Pu. 



TABLE VI: Correlation coefficients otij obtained from 2- 
dimensional Gaussian fits applied to the distributions from 
Figure [6] 





235 tj 


238-q 


239p u 


241 Pu 


'2'i'ojj 
238JJ 
239p u 

241 Pu 


1 

-0.29 
-0.62 
-0.48 


-0.29 
1 

0.38 
0.39 


-0.62 
0.38 
1 

0.84 


-0.48 
0.39 
0.84 
1 



all 159 fuel samples comparisons 
cients a 

2-dimcnsional Gaussians 



Correlation coeffi- 
y are obtained by fitting these distributions to 



Ae 



.i T -br(A?+A 3 2 -2a IJ A I A J ) 



(10) 



where Aj = (Sfi/fi—fii)/(Ti and A is a normalization con- 
stant. In the fit, fii and a% are allowed to float, resulting 
in fit values consistent in all cases with the values listed 
in Table [V] The fit values of otij are listed in Table IVII 

We find a weak anti-correlation between 235 U and all of 
the other isotopes. This is expected since 235 U dominates 
the thermal power production of the reactor, and so to 



IV. UNCERTAINTY IN THE ANTI-NEUTRINO 
SIGNAL CALCULATION 

We now examine how to propagate the uncertainties in 
Wth and the fi, found in the previous sections, into the 
anti-neutrino signal uncertainty. This is not straightfor- 
ward because of the significant correlation between SWth 
and the Sfi, as indicated by Equation[3J To first order, an 
increase in the thermal power by some factor implies an 
increase in the fa by the same factor. In fact, if we use the 
Sfi/ fi obtained from each spent fuel element discussed in 
the previous section to extract the thermal power uncer- 
tainty via standard error propagation on Equation [31 we 
obtain a distribution oi5W t h/W t h with width 2.5%, con- 
sistent with the expectation for a set of reactors instru- 
mented with Venturi-type flow meters of varying ages. 

The point of the thermal power normalization in Equa- 
tion [4] is to break this first-order correlation between the 
fi and Wth, and work with the fission fractions fi/F 
rather than the bare fi. But even the fi/F have a resid- 
ual correlation with the thermal power that is due to the 
fact that U isotopes are depleted while Pu isotopes are 
bred during reactor operation. An over-estimate of the 
power will lead to an over-estimate of the Pu production, 
and hence a higher estimated fraction of fissions from the 
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Pu isotopes. This is significant for the v e signal estima- 
tion because, by coincidence, the are slightly higher 
and the rii are slightly lower for the Pu isotopes than for 
the U isotopes. Since the et appear in the denominator 
of Equation [4] and the rii appear in the numerator, the 
effect is to lower slightly the signal estimate. 

Hence we see that the contributions to the signal un- 
certainty from the errors on Wth and the fi/F are slightly 
anti-correlated. It is this anti-correlation that is respon- 
sible for the slope of "Power" in Figure [2] being slightly 
less than one. However, the anti-correlation is only slight, 
and moreover implies a weak cancellation between the 
uncertainty contributions evaluated from the indepen- 
dent variation of Wth and the fi/F. To maintain sim- 
plicity in our analysis while remaining conservative in 
the error propagation, we ignore this anti-correlation and 
treat the errors due to these terms as being uncorrelated. 
We thus evaluate the uncertainty on the v e signal, a v , 
according to the simple quadratic sum of contributions 
from each of the terms appearing in Equation [3J 

°"n„ = °W + °/ + a 1 + °ither- ( 12 ) 

Since Equation [?] is linear in Wth, the first term, crw, is 
given by the thermal power uncertainty, SWth/Wth, dis- 
cussed in detail in Section [Til The contribution cr/ from 
the fission calculation uncertainties, Sfi/fi is more com- 
plex and will be described below. The third term, cr e , is 
the contribution due to the uncertainty in the e, . The last 
term, (J other, represents uncertainties in detector-specific 
components and other terms (such as the dN^/dE^) that 
have been explored elsewhere. This term will not be ad- 
dressed further here. 

Before addressing cr/, we briefly note that a e can be 
evaluated from error propagation on Equation^ Assum- 
ing that the errors on each of the are uncorrelated, we 
obtain 

The values of the and their uncertainties Sei are given 
in [H [65| as e( 235 U) = 201.7 ± 0.6 MeV, e( 238 U) = 
205.0 ± 0.9 MeV, e( 239 Pu) = 210.0 ± 0.9 MeV, and 
e( 241 Pu) = 212.4 ± 1.0 MeV. To evaluate <r e , we need 
to choose values for the fission fractions, fi/F. As exam- 
ples, Table I VIII shows the fission fractions for the three 
reactor units of PVNGS (PV1, PV2, and PV3), as well as 
averaged fission fractions of the reactors in Japan (KL), 
observed by the KamLAND experiment. The fission frac- 
tion values vary from reactor to reactor depending on 
the initial enrichment level and power history of the fuel. 
Table I VIII gives the corresponding value of a e , calculated 
with Eq. [T3l In all cases considered, the value of er e is 
-0.2%. 

Correlations between the Sfi / fi make an analytical es- 
timate of Of difficult. Instead we leverage the measured 
5fi/fi from Section UTTI to compute, for each of the 159 
spent fuel element comparisons, a value of cr/ according 



TABLE VII: Fission fractions i-, a e , and Of for the three 
PVNGS reactors (PV1, PV2, PV3), and for all reactors 
viewed by KamLAND 3|] (KL), averaged over the periods 
of operation. 



Reactor or 




h 

F 


■ [%] 




o c [%] 


°7 


[%] 


reactor group 


235jj 


238jj 


239p u 


241 Pu 




syst 


± stat 


PV1 


58.0 


7.4 


29.2 


5.4 


0.217 


0.13 


± 0.84 


PV2 


54.4 


7.5 


31.8 


6.3 


0.216 


0.13 


± 0.89 


PV3 


57.7 


7.4 


29.2 


5.7 


0.216 


0.13 


± 0.84 


KL 


56.3 


7.9 


30.1 


5.7 


0.216 


0.14 


± 0.88 



to 



i + (Ei£f "j)/(Ei£nj) 

(14) 

This equation may be obtained by replacing the fi in 
Equation^ with fi + Sfi, and n v with n v + Sn v . By sub- 
stituting the 159 values of Sfi/fi into Equation [Ml we 
maintain the implicit correlations between different iso- 
topes. And as explained above, the normalization by the 
thermal power removes the influence of the Wth uncer- 
tainty on the v e signal, so that cr/ accounts for the con- 
tribution from Sfi/fi only. To compute the values of the 
rii in Equation ll41 we assume no oscillation (i.e. short L) 
and an energy-independent efficiency above the inverse 
beta-decay threshold. For the fi/F we use the same sets 
of values used to evaluate cr e listed in Table IVIII 

An example cr/ distribution is drawn in Figure [7] us- 
ing the fission fractions from 0], listed in the "KL" row 
of Table IVTT1 When drawing the distribution of cr/, we 
weight the value from each spent fuel element compari- 
son according to its burnup to account for the fact that 
the higher-burnup samples near the center of the reac- 
tor contribute proportionally more to the v e signal than 
the lower-burnup samples near the core edges. The set 
of comparisons still includes a disproportionately large 
number of fuel elements from regions that are challeng- 
ing for the simulation; by including all such samples we 
ensure a conservative estimate of cr/. These problem- 
atic samples are responsible for the relatively large tails 
and the few outliers in the distribution. In order to in- 
corporate these samples without letting them dominate 
the characterization of the distribution, we fit the distri- 
bution to a single Gaussian, and take the fit mean and 
standard deviation as estimates of the systematic and 
statistical components, respectively, of 07. We repeated 
this procedure for each of the four choices of the fi/F 
listed in Table IVlTl As can be seen from these examples, 
cr/ is typically about 0.9%, of which ~0.1% is correlated 
between different reactors. 

For an experiment at a single reactor, unless the ex- 
perimenters and reactor operators have gone to extremes 
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H = 0.14 ±0.07 
a = 0.88 ± 0.05 




FIG. 7: o f distribution for the fuel comparisons, calculated 
with Eg. 1141 weighted by sample burnup and fit to a Gaussian. 
The mean of the fit is identified as the systematic component 
of (J f , while the spread represents the statistical uncertainty. 
The inclusion of a disproportionately large number of sam- 
ples from regions that are challenging to simulate contributes 
larger tails and a higher mean, leading to more conservative 
estimates uncertainties. 



to validate the core simulations, the uncertainty in the 
signal due to the error in the fission fraction calculations 
should be taken to be the full statistical uncertainty on <j/ 
added in quadrature with the systematic component. As 
an example, using the values listed in Table HT1 with Type 
I transit-time UFMs and incorporating the enthalpy un- 
certainties (as given in Table fl} to obtain the thermal 
power uncertainty, the resulting uncertainty in the anti- 



neutrino yield, ignoring a other, is 



'w 



'f 



« ^(0.2 2 -K).4 2 +0.15 2 +0.2 2 )+(0.14 2 +0.9 2 )+0.2 2 
« 1.1%. 

(15) 

If the reactor instead measures the thermal power output 
at the 2% level, then this uncertainty evaluates to 2.2%. 

For a multi-reactor experiment, some reduction of the 
statistical components of aw and Of is achieved. If re- 
actor k contributes a fraction p^ to the total reactor D e 
signal, then the total random error in the summed sig- 
nal is reduced by a factor of \/J2kPk- For example, 
consider the case of an experiment receiving an approxi- 
mately equal signal from 2 reactors equipped with Type 
I transit-time UFMs. The resulting uncertainty in the 
anti-neutrino yield (again, ignoring <J ther) is 



iL|l + 0.4 2 + 245l+0.2 2 + 0.14 2 + 1^|1 +0.2 2 



0.83%. 



TABLE VIII: Signal uncertainties at KamLAND, a v , ob- 
tained by combining aw, in quadrature with 07 and <r e . The 
calculations were performed assuming various power uncer- 
tainties. The 2.1% entry corresponds to the value used by 
the KamLAND collaboration. 



<" Elc [%] 


<TW [%] 


*,[%] 


*Y C [%] 


Te[%] 




systistat 


systistat 


systistat 






2.10±0.00 


2.10±0.00 


0.14±0.18 


1.0 


0.216 


2.3 


0.60±1.40 


0.60±0.26 


0.14±0.18 


0.2 


0.216 


0.76 


0.45±0.25 


0.45±0.05 


0.14±0.18 


0.2 


0.216 


0.59 



(16) 



As a more complex example, we study the thermal 
power and fission fraction uncertainties in the Kam- 
LAND experiment. We use the actual distance and nom- 
inal thermal power of the reactors to compute the frac- 
tions pk they contribute to the signal [23| . We then al- 
low all random uncertainties to be reduced by a factor 
°f vEii Pk ~ \ ■ F° r the fission fraction uncertainty, the 
0.88% statistical contribution for a single reactor is re- 
duced to 0.18%. However, for KamLAND we must add 
an additional uncertainty, crj^ , to account for the fact 
that to compute the fission fractions the KamLAND col- 
laboration used a generic reactor simulation which was 
found to agree with detailed simulations, such as those 
considered in this paper, at the 1% level (22J. In this 
reference, this uncertainty appears to be dominated by 
random errors, so it may be appropriate to allow this 
component to also be reduced by a factor of 

For the thermal power uncertainty in the KamLAND 
example, we consider three scenarios: the standard 
KamLAND assumption of an across-the-board 2.1% sys- 
tematic uncertainty, a slightly more realistic yet still 
conservative assumption that all Japanese reactors are 
equipped with properly-calibrated Venturi flow meters, 
and an ideal scenario in which all Japanese reactors are 
equipped with ultrasonic flow meters. The resulting val- 
ues of cr*° tal are listed in Table IVIIII For the standard 
KamLAND 2.1% thermal power error and with aj' IC 
taken to be the full 1%, the contributions of af and a e 
are negligible, giving a result consistent with the 2.3% 
uncertainty used by the KamLAND collaboration, and 
justifying a classification of the KamLAND treatment 
as "conservative" . However, if we make the assumption 
that Japanese reactors are equipped with flow meters 
at least as performant as Venturi flow meters, and fur- 
ther take the spread in crf IC evident in [12] to indicate 
that this term is predominantly random in nature, sig- 
nificant cancellations can be achieved. In this case, the 
systematic uncertainty due to the fouling is the dominant 
component; including all other components gives a total 
uncertainty of 0.76%. In an ideal scenario in which all 
Japanese reactors are equipped with Type I transit-time 
UFMs, the total uncertainty would drop to 0.59%. The 
dominant component is still the systematic uncertainty 
of the flow meter, but the relative contribution of the 
other terms is significant. 

The actual situation in the Japanese reactor industry is 
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probably somewhere in between the all-Venturi case and 
the ideal case of all UFMs. Many reactors may already 
be equipped with UFMs, but some may not, and those 
that are may simply use them to calibrate their Venturi 
meters. However, the total uncertainty, a nv , does not 
change much between these two cases. Thus for Kam- 
LAND, a reduced error of 0.6-0.8% may be more ap- 
propriate. A better understanding of the nature of the 
Japanese instrumentation and reactor simulations may 
result in yet smaller overall uncertainties for KamLAND. 

V. CONCLUSION 

The uncertainty in the estimated anti-neutrino inter- 
action rate at reactor experiments has contributions from 
the thermal power and the fission calculation uncertain- 
ties. The dominant contributions to the thermal power 
uncertainty come from the enthalpy rise and mass flow 
rate in the steam generator (PWR) or reactor vessel 
(BWR). The former has a systematic uncertainty of 0.2% 
and a random uncertainty of typically 0.15%. The uncer- 
tainty in the mass flow rate depends on the instrumenta- 
tion used to perform the measurement. Traditional Ven- 
turi flow meters exhibit significant fluctuations of 1.4% 
and a systematic uncertainty of 0.6% due to fouling. 
More precise ultrasonic flow meters have been evaluated 
to have systematic uncertainties below 0.5% and typically 
show random fluctuations below this level. 

We estimated, for the first time to our knowledge, the 
contribution to the v e signal uncertainty from errors in 
the fission rate calculations. These errors were extracted 
from comparisons of heavy element concentration mea- 
surements in spent fuel elements with calculations of 
those concentrations by the same simulations that re- 
actor v e experiments rely on to obtain the fraction of 
fissions from each of four main heavy isotopes. Error 



propagation accounting for correlations between the four 
isotopes yielded a 0.1% systematic uncertainty on the v e 
signal rate, and a 0.9% random uncertainty. We also es- 
timated a 0.2% contribution due to the uncertainty in 
the measured energy release per fission of each isotope. 

We demonstrated how these errors can be combined 
for typical reactor experiment configurations, and high- 
lighted situations in which cancellations occur between 
different reactor sites. With this methodology, we found 
that these contributions to the signal rate uncertainty in 
the KamLAND experiment may be reduced from 2.3% 
to 0.76%, or even 0.59%, pending better understand- 
ing of reactor instrumentation. As we move deeper into 
the phase of precision reactor neutrino oscillation exper- 
iments, and in particular multi-reactor #13 experiments, 
such sub-percent understanding of these errors and their 
systematic and random components will become more 
and more important. As reactor engineering, flow rate 
metrology, and core modeling continue to improve, fur- 
ther reductions may be achievable in these uncertainties. 
Toward this end, closer relationships between neutrino 
scientists and reactor operators should be encouraged. 
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